Root-based inorganic carbon uptake increases the growth of Arabidopsis thaliana and changes transporter expression and nitrogen and sulfur metabolism

Root-based uptake of inorganic carbon has been suggested as an additional carbon source. Our study aimed to characterize and understand the root-based uptake and fixation mechanisms and their impact on plant growth. 13C-labeled bicarbonate fed to Arabidopsis roots was assimilated into aspartic acid but mainly into sucrose, indicating that the added inorganic carbon was transported to the leaves. A hydroponic treatment was also established for A. thaliana using 2 mM NaHCO3 at pH 5.6, which enhanced the photosynthetic and growth parameters. According to transcriptome sequencing data, the observed enhancement in growth may be orchestrated by trehalose-6-phosphate signaling and supported by augmented nitrogen and sulfur assimilation. The analysis also revealed regulatory and transporter activities, including several nitrate (NRT2.1), and sulfate transporter (SULTR1;1 and SULTR1;2) candidates that could participate in bicarbonate uptake. Different transporters and carbon fixation mutants were assessed. Arabidopsis homologs of SLOW-TYPE ANION CHANNEL 1 (slah3) CARBONIC ANHYDRASE (βca4), and SULFATE TRANSPORTER (sultr1;2) mutants were shown to be inferior to the bicarbonate-treated wild types in several growth and root ultrastructural parameters. Besides, aquaporin genes PIP1;3 and PIP2;6 could play a negative role in the carbon uptake by venting carbon dioxide out of the plant. The findings support the hypothesis that the inorganic carbon is taken up by the root anion channels, mostly transported up to the shoots by the xylem, and fixed there by RuBisCo after the conversion to CO2 by carbonic anhydrases. The process boosts photosynthesis and growth by providing an extra carbon supply.


Introduction
Increasing crop yields is a crucial breeding objective, considering global climate change, overpopulation, and the use of crop plants as energy sources.Various studies have suggested that one of the most practical ways to boost crop yield is to increase photosynthetic activity (E ́va et al., 2019;Furbank et al., 2020).Promising new research aims to make RuBisCo more specific for CO 2 (carboxylation reaction) and to increase the speed of the enzyme by improving its kinetics (Gunn et al., 2020;Lin et al., 2022;Zhou et al., 2023).Moreover, several studies have explored ways to enhance CO 2 supply to RuBisCo in C 3 plants, such as implementing C 4 photosynthesis into rice (C 4 Rice Project) (Von Caemmerer et al., 2012;Karki et al., 2013;Furbank et al., 2023), introducing a photorespiratory bypass within the chloroplast of Arabidopsis (Kebeish et al., 2007), tobacco (South et al., 2019), and rice (Shen et al., 2019), and ectopically over-expressing algal bicarbonate transporters in tobacco (Nölke et al., 2019).Considerable work has been done to better understand leaf CO 2 uptake and explore possibilities for photosynthetic improvement.However, it has long been known that the rhizosphere contains more CO 2 than the atmosphere (Nel and Cramer, 2019).Currently, the global average concentration of atmospheric CO 2 as of June 2024 is 426.91 ppm (http://www.co2.earth/), while Sotomayor and Rice (1999) reported the annual mean soil CO 2 concentration at a depth of 3 m to be 11000 ppm, Yonemura et al. (2009) reported 1000-23000 ppm at a depth of 10 cm with higher values measured in non-tilled soils.Even at equal levels, CO 2 uptake in the damp environment of the roots could possibly provide the benefit of less water loss for the plant compared to the leaves with open stomata.
E ́va et al. (2019) suggested that root-based inorganic carbon uptake and subsequent transport to shoots would increase photosynthetic fixation.Although the exact mechanism remains unclear, this process has been observed in plants.Different examples of carbon fixation, such as the transport and recycling of root-respired carbon, have been proposed as a means of coping with drought stress (Søndergaard and Sand-Jensen, 1979;Vapaavuori and Pelkonen, 1985;Biatczyk et al., 1994;Bloemen et al., 2013;Dabrowska-Bronk et al., 2016;Han et al., 2022).Moreover, Dabrowska-Bronk et al. (2016) confirmed that soluble bicarbonate ions were absorbed from the soil solution by the root hair.As well, it has been demonstrated that increased sucrose concentrations in plants exposed to high CO 2 levels affect root development.This and other evidence showing ultrastructural changes in the roots of high-carbon-treated plants indicate that root anatomy could play a role in carbon uptake.Similarly, it has been observed that high CO 2 levels increase the total number, diameter, and length of Arabidopsis thaliana roots (Thompson et al., 2017).Furthermore, it has been demonstrated that using lower bicarbonate concentrations 0.7 mM (Vapaavuori and Pelkonen, 1985), 1 mM (Wanek and Popp, 2000), 2 mM (Li et al., 2023), 3 mM (Dabrowska-Bronk et al., 2016), and 5 mM (Biatczyk et al., 1994) could promote growth, improve intracellular water metabolism, nutrient transport, and photosynthetic capacities.Conversely, higher concentrations such as 6 mM (White and Robson, 1990), 7 and 12 mM (Li et al., 2023), 10 mM (Zribi and Gharsalli, 2002), and 20 mM (Alhendawi et al., 1997) can inhibit plant growth, impair iron acquisition, inhibits effects on plants photosynthetic rate, stomatal conductance, and transpiration rate.
The pH of the solution determines the chemical state of inorganic carbon.At pH 5, most of the inorganic carbon is found in the form of CO 2 , whereas at pH 8, the majority is in the form of bicarbonate ions (Shimono et al., 2019).There are conflicting results regarding the effects of high and low pH values, and further research is required to determine the optimal pH for inorganic carbon uptake.Some studies have shown increased absorption and root-to-shoot transit of inorganic carbon in beans and sunflowers at a low pH (Shimono et al., 2019).For instance, when A. thaliana was irrigated with a bicarbonate solution at pH 7 that is, a high pH solution rich in bicarbonate-it increased plant development (Dabrowska-Bronk et al., 2016).
Energy-driven uptake of bicarbonate ions is a potential mechanism for root-based inorganic carbon uptake.To the best of our knowledge, true bicarbonate pumps have not yet been isolated from land plants.However, bicarbonate ions were found to inhibit nitrate, sulfate, and phosphate uptake, suggesting a competitive mechanism and indicating that anion transporters with low specificity might also take up bicarbonate (Poschenrieder et al., 2018).Besides the anion transporters, the potential involvement of CO 2 -permeable aquaporins in carbon uptake cannot be ignored, especially at low pH (CO 2 -rich environment).The high protein and sterol content of plant membranes typically act as barriers to CO 2 diffusion (Groszmann et al., 2017).CO 2permeable aquaporins, such as Arabidopsis AtPIP1;2 and barley HvPIP2;1, are highly expressed in the roots (Katsuhara et al., 2002;Postaire et al., 2010).Their function was previously shown to assist CO 2 in freely crossing membranes and reaching chloroplasts in leaves (Uehleln et al., 2003;Hanba et al., 2004;Heckwolf et al., 2011;Ermakova et al., 2021).Although aquaporins have been shown to play a role in water uptake, their potential contribution to rootbased CO 2 acquisition is yet to be studied (Katsuhara et al., 2002;Postaire et al., 2010).
Following root-based carbon uptake, the plant may fix the carbon or release it back into the environment (Shimono et al., 2019).Carbonic anhydrases (CA) catalyze the reversible conversion between the inorganic carbon forms, which can then serve as a substrate of RuBisCo (CO 2 ) or PEP-carboxylase (HCO 3 -).All three types of CAs (a-, b-, and g) are present in plants in various organs (leaf, seed, root) at various subcellular locations (chloroplast, mitochondrion, cytoplasm, plasma membrane), with the b-CAs being the most abundant (DiMario et al., 2017).Traditionally it was thought that the main role of chloroplastic CAs in mesophyll cells of C 3 plant leaves is to supply RuBisCo with CO 2 but it found only weak experimental support, except for seedlings (Ferreira et al., 2008).However, CAs in leaves of C 3 plants take part in various processes like regulation of stomatal movement, biotic and abiotic stress responses, amino acid biosynthesis, and lipid biosynthesis (DiMario et al., 2017;Rudenko et al., 2022).Of the Arabidopsis bCAs, AtbCA4.2(located in the cytoplasm) and AtbCA5 (located in the chloroplast) showed the highest expression in roots (Dimario et al., 2016).A longer isoform (alternative transcript) of AtbCA4 (AtbCA4.1) is specifically expressed in the plasma membrane of leaf cells, and together with the apoplastic aCA2 plays an important part in the transfer of CO 2 and HCO 3 − to the plant cell (Dimario et al., 2016;Weerasooriya et al., 2024).PEPC catalyzes the carboxylation of PEP in the presence of HCO 3 -and Mg 2+ to yield oxaloacetate and inorganic phosphate.In C 4 and CAM plants, the photosynthetic tissues contain high levels of PEPC that catalyzes the initial fixation of atmospheric CO 2 during photosynthesis.The much lower levels of PEPC seen in the leaves of C 3 plants contribute to the enzyme's anaplerotic function with a key role in the coordination of C and N metabolism, replenishing the TCA cycle with intermediates used for amino acid biosynthesis, and play a role in the regulation of the cellular pH (Miyao and Fukayama, 2003;Nimmo, 2003).It has been described that the fixation of DIC (Dissolved Inorganic Carbon) provided for the roots is associated with the activity of the PEPC (Bialczyk and Lechowski, 1995).Some examples of the PEPC localized in the roots fixing CO 2 are described byHibberd and Quick (2002) in tobacco and by Brown et al. (2010) in Arabidopsis.
Evidence from the literature supports root-based inorganic carbon uptake, as was previously mentioned.The process behind carbon uptake and fixation by the roots, as well as how it affects plant growth and homeostasis, is still poorly understood.It is also unclear whether carbon is predominantly fixed in the roots and exported to the shoots in an organic form or transported in an inorganic form and fixed in the leaves.Our long-term hypothesis is that our crop plants' photosynthetic performance and osmotic stress tolerance can be improved by enhancing root-based inorganic carbon uptake.The present study is the first, screening part of the project, aiming to establish a growth-promoting inorganic carbon treatment for Arabidopsis and use it to characterize the root-based uptake, fixation, and growth promotion in terms of anatomy, physiology, and molecular biology.A better understanding of these processes during the present study shall lead to the identification of candidate genes like possible bicarbonate transporter candidates which can be validated and tested in future works and over-expressed in transgenic plants, aiming to boost their photosynthesis and osmotic stress tolerance.

Growth tests
Arabidopsis seeds were surface-sterilized first with 1 ml of 70% (v/v) ethanol (freshly diluted from absolute ethanol supplied by VWR International, Radnor, PA, USA) and then by rotating for 12 min in a mixture of 50% (v/v) ethanol, 1.5% bleach (Feel-It, Egyhaźasgerge, Hungary), and 0.05% (v/v) Tween 20 (Merck-Sigma Group, Darmstadt, Germany).The seeds were subsequently rinsed with 96% (v/v) ethanol and washed four times with sterile water.Sterile seeds were cold-stratified (4°C) for two days in the dark.The seeds were then germinated, and the plants were grown in solid ½ MS medium (Murashige and Skoog Basal Salt Mixture, pH 5.7; Merck-Sigma Group) with 1% (w/v) sucrose (VWR International) in vertically kept square Petri dishes until they had a well-developed root system (2 weeks) Supplementary Figure 1A.In the third week, the treatment started for the standardization of the growth promotion system (Supplementary Figure 1B).The plants were transferred to hydroponic cultures supplied with ¼ MS salts together with two different pH: 5.6 and 7.2 and three different concentrations of NaHCO 3 (Merck-Sigma Group) (0, 2, 4 mM).Once the growth promotion was standardized the mutant's evaluation came (Supplementary Figure 1C), the treatment used was ¼ MS salts at pH=5.6 and 2 mM of NaHCO 3 , which was the system with the best results in the standardization process.The hydroponic system was similar to that described by Tocquin et al. (2003) but with the addition of an aquarium air pump.The hydroponic solution was changed weekly as mentioned in Robison et al. (2006); Conn et al. (2013).

Physiological measurements
General growth indicators, such as root length (RL) and rosette diameter (RD), were measured using a ruler, and total fresh root weight (RW) and total fresh rosette weight (ROW) were measured on a scale Sartorius CP4202S.

Gas exchange analysis
Leaf gas exchange analysis was performed on the third fully developed leaf of 4-week-old plants one day after the last treatment with NaHCO 3 .The CIRAS-3 portable photosynthesis system (PP Systems, Amesbury, MA, USA) was used with a narrow leaf chamber (1.75 cm 2 ).CO 2 fixation rate of control and 2 mM NaHCO 3 treated plants were investigated at 390 ppm CO 2 in an infrared gas analysis cuvette, and by maintaining leaf and cuvette temperatures at 21°C.The actinic light was provided by an LED unit of CIRAS 3 (PPFD = 250 µmol m -2 s -1 , RGBW = 60:5:25:10).The net photosynthetic rate (P n , µmol CO 2 m -2 s -1 ), stomatal conductance (gs, mmol CO 2 m -2 s -1 ), transpiration rate (E, mmol H 2 O m -2 s -1 ), and intracellular CO 2 concentration (Ci, µmol CO 2 mol air -1 ) were recorded at the steady-state level of photosynthesis (5-min-lag-time was allowed, and data were recorded at the 15 th min of the measurement).

Fluorescence imaging analysis
Chlorophyll-a fluorescence induction (FI) analysis was done in the plants grown during the standardization of the growthpromoting system (0, 2, 4 mM NaHCO 3 and two different pH 5.6 and 7.2).The analysis was carried out using a pulse amplitude modulated fluorometer (PAM) with a blue LED-Array Illumination Unit IMAG-MAX/L (l=450 nm) (Imaging-PAM MSeries, Walz, Effeltrich, Germany).After 30-min dark adaptation of the plants, the F 0 , F m, and F v /F m parameters were determined using a saturation pulse (SP, PPFD = 3000 µmol m -2 s -1 with 0.8 s duration).The quenching analysis was carried out after a 40-sec lag time at laboratory temperature (adapted state) using continuous blue actinic light (PPFD = 250 µmol m -2 s -1 ) and 30-sec SP frequency until the steady-state level of photosynthesis was reached (total duration: 15 min).The actual PSII quantum yield [Y(II)], the quantum yield of the regulated way of energy dissipation [Y(NPQ)], and the quantum yield of the non-regulated way of energy dissipation [Y(NO)] parameters were calculated during analysis, as described by Klughammer and Schreiber (2008).

Fluxomic analysis
For the fluxomic analysis xylem sap and phloem sap were collected from Col plants grown in hydroponic cultures with the standardized growth promotion system (control; 0 mM NaHCO 3 and treatment; 2 mM NaHCO 3 and pH 5.6).The plants used were four old weeks and the administration of the labeled carbon ( 13 C NaHCO 3, Merck-Sigma Group) was added one day before the sample collection.The treatment included three groups: control (not receiving any NaHCO 3 ), 12 C-treated (receiving 2 mM 12 C NaHCO 3 , purchased from Merck-Sigma Group, carbon composition according to natural abundance: 1.1% 13 C), and 13 Ctreated (receiving 2 mM 13 C NaHCO3, purchased from Merck-Sigma Group, containing 99% 13 C).In the case of the collection of xylem sap, we followed the methodology of Li et al. (2010); and for the phloem sap, we followed the methodology of Guelette et al. (2012).
From the phloem samples, 300 µL aliquots were evaporated in a vacuum centrifuge and then reconstituted in 300 µL of 9:1 (v/v) acetonitrile: water solution containing 0.5% (v/v) formic acid.The reconstituted samples were vortexed and then centrifuged at 16,000 g for 15 min at 4 °C; afterward, 125 µL aliquots were taken from the supernatants and filled into LC-MS vials equipped with microinserts.From the xylem samples, 15 µL aliquots were diluted with 135 µL of acetonitrile containing 0.55% (v/v) formic acid.The samples were vortexed and then centrifuged at 16,000 g for 15 min at 4 °C; afterward, 100 µL aliquots were taken and handled as above.
To determine the multiple reaction monitoring (MRM) transitions for the quantification of labeled and non-labeled analyses, ESI-MS/MS analysis was carried out using an Acquity Iclass ultra-performance liquid chromatography (UPLC) system coupled to a Xevo TQ-XS Triple Quadrupole Mass Spectrometer (Waters, Milford, MA, USA), used in electrospray ionization (ESI) negative mode.Separation was performed on a ZIC-cHILIC column (Merck; 100 mm * 2.1 mm * 3.0 µm) at 40 °C.For gradient elution, water and acetonitrile containing 0.5% (v/v) formic acid were used.The applied gradient and the UPLC-ESI-MS/ MS parameters (Supplementary Table 2) are presented together with the relevant MS/MS spectra (Supplementary Figures 2-6).For the selective quantification of non-labeled and [1-13 C glc ]-labeled sucrose molecules, the MRM transitions of m/z 341→161 and m/z 343→163 were used, respectively.The intensity values of m/z 343→163 ('A+1' isotopolog for [1-13 C glc ] sucrose in ESI negative mode) determined for the samples in the labeling experiment were corrected with the residual m/z 343→163 data of non-labeled sucrose standard ('A+2' isotopolog for natural sucrose in ESI negative mode).
For the selective quantification of natural aspartic acid and aspartic acid-4-13 C molecules, the MRM transitions of m/z 132→88 and m/z 133→88 were used, respectively.The intensity values of m/z 133→88 ('A' isotopolog for aspartic acid-4-13 C in ESI negative mode) determined for the samples in the labeling experiment were corrected with the residual m/z 133→88 data of non-labeled aspartic acid standard ('A+1' isotopolog for natural aspartic acid in ESI negative mode).

Light microscopy
The required chemicals were from Merck-Sigma Group, unless otherwise stated.For root anatomy, the following genotypes were evaluated Col, pepc3, bca4, sultr1;2, slah3, pip1;2-1, pip1;2-2 and the treatment was 0 mM or 2 mM NaHCO 3 the samples were obtained from the root differentiation zone with root hairs.The roots were fixed in 50 mM Na-cacodylate buffer (pH 7.2) containing 2.5% (v/v) glutaraldehyde, and 4% (v/v) formaldehyde overnight at 4°C, washed, dehydrated in an ethanol series, and gradually infiltrated with London Resin (LR) White acrylic resin (Ted Pella, Redding, CA, USA) according to the manufacturer's instructions.The resin was polymerized under UV light at -20°C.Semi-thin sections (1 mm) were serially sectioned at the transverse plane of the roots using an Ultracut-E microtome (Reichert-Jung, Heidelberg, Germany) and stained with periodic acid-Schiff (PAS) and 1% (w/v) Amido Black for polysaccharides and proteins, respectively.The stained sections were mounted in 50% (v/v) glycerol containing 7% (v/v) acetic acid and examined under a BX51 light microscope (Olympus, Tokyo, Japan).Photographs were captured using a TrueChrome-II color camera (Tucsen Photonics, Gaishan Town, Fujian, PR China) with integrated software to measure the parameters of root diameter (RTD), epidermis width (EpWd), endodermis width (EnWd), cortex width (CxWd), stele diameter (StDm) and outer cortical cell width (OuCCWd).The anatomy and various parts of the Arabidopsis roots were identified according to Buchner et al. (2004).

Transcriptome sequencing
Samples from middle-aged rosette leaves and roots of approximately 4 weeks old plants were collected at the end of the bicarbonate treatment.For total RNA extraction, 70 mg of each sample in three repetitions were homogenized in 500 mL of TRI-Reagent (Zymo Research, Irvine, CA, USA).Total RNA was extracted using the Direct-zol RNA MiniPrep System (Zymo Research), according to the manufacturer's protocol, including the on-column DNase digestion.IRNA concentrations were measured using a NanoDrop 2000 Spectrophotometer (Thermo Scientific), and the A260/A280 ratio was used to estimate the RNA quality.Total RNA samples were then subjected to quality control and custom sequencing.All our samples passed the quality control and had an RNA integrity number (RIN) of ≥7.The general parameters of the analysis included the production of 150-bp pair-end reads and 12 Gbp raw data per sample on the Illumina NovaSeq 6000 platform (Novogene, Nanjing, China), and the reads were mapped to the ensemblplants_arabidopsis_thaliana_ tair10_gca_000001735_1 genome assembly.Only differentially expressed genes (DEGs) with significant (p<0.05) and log 2 FC ≥ | 1| were involved in any further investigation, including gene ontology (GO) enrichment (Ashburner et al., 2000) and Kyoto Encyclopedia of Genes and Genomes [KEGG; (Kanehisa and Goto, 2000)] pathway (org code: ath) analysis.These functional analyses have been implemented by the clusterProfiler R package, in which gene length bias was corrected.The graphics for the volcano plots were done using Venny 2.1 (https://bioinfogp.cnb.csic.es/tools/venny/index.html),the heatmaps were prepared with Heatmapper (http://heatmapper.ca/),and the bar graphics and the enrichment bubble plot with SRplot (https://www.bioinformatics.com.cn/en).Four treatments are reflected in the analysis, the MSRC (Murashige and Skoog Root Control), MSRT (Murashige and Skoog Root Treated with the 2mM NaHCO 3 ), MSLC (Murashige and Skoog Leaf Control), MSLT (Murashige and Skoog Leaf Treated with the 2mM NaHCO 3 ).

qPCR validation of the RNAseq data
Total RNA was extracted in the same way as described above.The samples (1,000 ng in each reaction) were then reversetranscribed using random primers and the RevertAid First Strand cDNA Synthesis Kit (Thermo Scientific) according to the manufacturer's instructions.The obtained cDNA was diluted threefold and used as a template for qPCRs.The reactions were run in triplicate using the CFX96 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA) and the PCRBIO SyGreen Mix (PCR Biosystems, London, UK).
Several transporter candidates and genes participating in cell wall composition were amplified, and their expression was compared between qPCRs and transcriptome sequencing.The transcription levels of the studied genes were calculated using the 2 -DDCt method (Livak and Schmittgen, 2001) with an efficiency correction step applied according to (Pfaffl, 2004).Primers used for qPCR validation are listed in Supplementary Table 3.

Statistical analysis
For the optimization of the growth promotion system, sample number (no. of plants measured) for the growth parameters was n=40, for the gas exchange analysis n=4, and for the fluorescence imaging n=10.This was individual for each of the group treatments of the pH: 5.6 and pH 7.2, with the respective treatments 0mM (control), 2mM and 4mM NaHCO 3 .The data is shown as mean ± SD and were statistically analyzed using a two-way ANOVA testing the effect of pH and the concentration of bicarbonate on the given growth or photosynthetic parameter followed by Tukey's HSD post-hoc test in RStudio 4.3.0.
For the evaluation of the mutant plants, sample number (no. of plants measured) for the growth parameters was n=12, for the gas exchange n=4, for the fluorescence imaging n=12 and for the light microscopy evaluation n=12.This is for each of the mutants evaluated with the corresponding treatment (0mM and 2mM).An independent sample Student's t-test (p<0.05) was performed for the comparison of the treatments per genotype in Microsoft Excel.Data are shown as the mean ± SD.
For the fluxomics analysis, six samples were evaluated using an independent sample Student's t-test (p<0.05) to compare the control and the treatment with labeled and non-labeled 2 mM NaHCO 3 and pH 5.6.

Optimization of growth-promoting inorganic carbon treatment for A. thaliana
Evidence from the literature suggests that low-concentration bicarbonate solutions have growth-promoting effects.Our experiments began with the growth optimization of hydroponic treatment in Arabidopsis thaliana (Col).Low (2 mM) and high (4 mM) NaHCO 3 treatments were tested at different pHs (5.6 and 7.2).We observed growth promotion (compared to untreated plants) at a 2 mM pH 5.6 NaHCO 3 condition.It improved the root and rosette weight and rosette diameter.Additionally, the net photosynthesis rate increased significantly at 2 mM NaHCO 3 and pH 5.6 (see Supplementary Figures 7-9).On the other hand, high-concentration (4 mM) bicarbonate treatments were inhibitory, especially at high pH to the point that gas exchange parameters could not be determined due to the state of the plants.

Uptake and fate of inorganic carbon supplied to roots
A fluxomic study was performed using carbon isotope ( 12 C and 13 C)-labeled NaHCO 3 to monitor the uptake of inorganic carbon and fixation of the derivative carbon in wild-type plants grown in the optimized, growth-promoting pH 5.6 hydroponic solution containing 2 mM NaHCO 3 .Aspartic acid, which can be produced by the transamination of oxaloacetate catalyzed by PEPC, was used as a marker for PEPC-based carbon fixation.Whereas sucrose is the main photoassimilate in land plants and can be considered an indirect marker of RuBisCo-based carbon fixation.The analysis revealed that 13 C-labeled sucrose could not be detected in the xylem sap (Figure 1A).However, 0.06 nM of 13 C of aspartic acid was found in the xylem sap (Figure 1B) of plants treated with 13 C, a three-fold increase compared to the control.In contrast to the xylem sap, but in line with the physiological function of the phloem, a considerable amount of 13 Clabeled sucrose was detected in the phloem sap (1.02 nM, a 20-fold increase in the concentration as a result of the 13 C treatment) (Figure 1A), suggesting that NaHCO 3 was taken up by the roots and the carbon was assimilated and incorporated into sucrose.The treatment also increased the ratio of 13 C sucrose compared to the 12 C one by 12-fold, reaching 2% of all sucrose (Supplementary Figure 10).Similarly, a 4-fold higher level of 13 C-labeled aspartic acid (0.22 nM) was detected in the phloem sap of treated plants, (Figure 1B), but its ratio reached only 0.4% of all aspartic acids.
Therefore, we proved the uptake and fixation of the hydroponically supplied carbon in Col during the optimized 2 mM NaHCO 3 treatment at pH 5.6, and we mostly found fixed carbon in sucrose in the phloem sap.

Transcriptional responses to NaHCO 3 treatment in roots and leaves in Arabidopsis plants
Although numerous transporters and channels are thought to be associated with the uptake and translocation of inorganic carbon, it is still unclear whether they are truly involved in CO 2 /HCO 3 uptake by roots and transport in plants.To identify more candidate genes involved in CO 2 /HCO 3 -translocation and to better understand the mechanism of fixation and the effect of carbon uptake, we performed RNA-Seq analysis of the roots and rosette leaves of 4-week-old wild-type Col-0 plants grown under control and 2 mM NaHCO 3 regimens in triplicate.B a s e d o n l o g 2 F C≥ 1 (upregulation), Log 2 F C ≤ -1 (downregulation) and p-adj ≤ 0.05, we identified 341 genes that were differentially expressed in the leaves of the control and treated plants.The number of differentially expressed genes (DEGs) that were upregulated in the leaves was 236 and downregulated by 105 (Figures 2A, C).In the roots, 1472 genes were differentially expressed after bicarbonate treatment.The number of upregulated DEGs was 615 and 857 were downregulated (Figures 2B, C).Overall, more genes were affected in the roots than in leaves.
Venn diagrams were used to investigate the presence of shared DEGs between the roots and leaves of control and treated plants.When comparing the DEGs reported for MSRC vs. MSRT and MSLC vs. MSLT (Figure 3A), only 80 DEGs were present in both roots and leaves.Furthermore, the results indicated that only 23 and 41 DEGs were commonly upregulated and downregulated, respectively, in the roots and leaves (Figures 3B, C).Additionally, 261 DEGs were exclusively present in MSLC vs. MSLT, with 82 upregulated and 195 downregulated genes (Figure 3D).For MSRC vs. MSRT, 1392 DEGs were expressed, with 834 upregulated and 574 downregulated genes (Figure 3D).A heatmap representation of To further investigate the function of the DEGs regulated in MSLC vs. MSLT and MSRC vs. MSRT, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were carried out.According to the GO analysis, 80 GO terms were assigned to the MSLC vs. MSLT group, while a total of 237 GO terms were significantly observed in the MSRC vs. MSRT group.There were 37 GO terms assigned to MSLC vs. MSLT and MSRC vs. MSRT.Some of the common GO terms were response to nitrogen compound (GO:1901698), trehalose biosynthetic process (GO:0005992), and response to starvation (GO:0042594).From upregulated DEGs, 49 GO terms were assigned to the leaves and 92 GO terms were assigned to the roots, of which three GO terms were common between the two: response to nitrogen compound (GO:1901698), response to reactive oxygen species (GO:0000302), and response to oxidative stress (GO:0006979).For the downregulated DEGs, 48 GO terms were assigned to the leaves, and 168 GO terms were assigned to the roots.Among the downregulated GO terms, nine were common to the leaves and roots, notably e.g.iron ion transport (GO:0006826) and cutin biosynthetic process (GO:0010143), (GO terms can be found in Supplementary Tables 4-7). Figure 4 shows the GO terms with a higher number of genes.
For the roots, the upregulated GO terms (Figure 4C) for BP were chemical homeostasis (GO:0048878), metal ion transport (GO:0030001), response to nitrogen compound (GO:1901698), and inorganic ion transmembrane transport (GO:0098660).For MF cation transmembrane activity (GO:0008324), and inorganic cation transmembrane transporter activity (GO:0022890).Interestingly, many GO terms related to transmembrane transporter activity, homeostasis, and cellular responses to stimuli were upregulated in the roots.The downregulated GO terms for the roots (Figure 4D) of BP were secondary metabolic process (GO:0019748), response to toxic substance (GO:0009636), and response to oxidative stress (GO:0006979).In contrast to the leaves, the presence of the GO term for the apoplast (GO:0048046) was downregulated.MF terms related to tetrapyrrole binding (GO:0046906) and heme binding (GO:0020037) contained the highest downregulated number of DEGs.
Most DEGs upregulated in the roots were related to transmembrane transport activity and homeostasis.We then decided to look forward to the presence of all transporters and channels.Tables 1, 2 list the up-regulated and down-regulated transporters and channels, respectively.We observed that sugar, ABC, vacuolar iron, and phosphate transporters were the families  with the most up-regulated DEGs.The nitrate transporter NRT2.1, together with the zinc transporter ZIP10, showed the highest foldchange; however, the expression of ZIP10 was not high.Similar to this upregulation, some ABC and sugar transporters were downregulated.The reliability of transcriptome sequencing was confirmed using qRT-PCR, which investigated several DEGs, including those with suspected involvement in root-based carbon uptake and fixation (NRT2.1,PEPC3, SULTR1;1; SULTR1;2; SLAH3, bCA4) between control and treated roots, as well as between leaf samples.As shown in Supplementary Figure 12, the qRT-PCR results were largely consistent with the transcriptome data.Furthermore, linear regression analysis was conducted to assess the correlation between qRT-PCR and RNA-seq, yielding R 2 values of 0.9564 for the roots and 0.9164 for the leaves, indicating a strong positive correlation between the two.

Physiological evaluation of possible inorganic transporters
After optimizing a growth-promoting treatment and proving the uptake and fixation of the added carbon, as well as identifying gene candidates that could take part in the uptake and fixation of the carbon, their gene knockout mutants were assessed under the growth-promoting system (2 mM NaHCO 3 at pH 5.6) to validate their function.
The data analyzed for the screening of candidates to be evaluated were based on a p-value of 0.05 and a log 2 FC of 0, genes with high expression levels, even if their FC did not change less than two-fold between control and treated plants, were studied based on literature information.DEG-lists are available as Supplementary Tables 12 and 13, our candidate list is as Table 3.
Growth parameters, photosynthetic parameters, and root microscopy results were evaluated for both the wild-type and mutant plants.All growth parameters including rosette diameter (ROD) (Figure 6A), rosette weight (ROW) (Figure 6B), root length (RTL) (Figure 7A), and root weight (RTW) (Figure 7B) were significantly increased in treated Col-0 plants compared to the control, whereas distinct trends were observed in the mutant lines.
Regarding the shoot system, the ROD of Col-0 exhibited a substantial growth of 27%.Mutants such as pepc3 and aquaporin displayed a notable increase, with aquaporins nearly doubling the growth percentage by 50% between the control and the treatment.On the other hand, bca4 showed a significant decrease in ROD, (approximately 9%).The slah3 and sultr1;2 mutants did not exhibit any significant differences as a result of the bicarbonate treatment (Figure 6A).The ROW of Col-0 demonstrated a substantial increase of 57%.The aquaporin mutants exhibited a significant increase of over 100%, particularly pip1;2-2 with more than a 500% difference, and slah3 also showed a significant increase of 60% (Figure 6B, Supplementary Table 14).Apart from comparing the treated plants in each genotype to their own control, it was also worth comparing the treated plants of different genotypes, here we could see that pip1;3, and pip2;6 had a bigger ROD and ROW compared to treated Col-0.
The shoot system of the plants were not just compared by measuring growth parameters, photosynthetic analysis was also conducted.We measured gas exchange to better understand the impact of the 2 mM NaHCO 3 treatment.Col-0 showed an increase in all the parameters assessed, except for the Ci, which decreased by 9% (Figure 6E).P n was significantly increased in the treated bca4, sultr1;2, and aquaporin mutants, except for pip1;2-1 (Figure 6C).While P n in sultr1;2 was significantly higher than the control, with a difference of 57%, the other evaluated parameters decreased compared with the control, but the difference was not statistically significant.Similar results were obtained for slah3 and pip2;6.Notably, bca4 was the only mutant that exhibited a significant increase in E and gs of 17% and 25%, respectively, similar to the wild-type (Figures 6F, D) (Supplementary Table 15).In a  6G, I) and treated (Figures 6H, J) Col-0 and the mutant lines are provided.
In terms of root measurements, the wild-type exhibited a significant 15% increase in RTL compared to the control, while pepc3, and sultr1;2 showed significant decreases of 13% and 11%, respectively (Figure 7A).Regarding RTW, Col-0 showed a significant increase of 88% with the treatment, except for pepc3, where all the mutants increased their RTW.The aquaporin mutants displayed a significant increase of over 140% (Figure 7B) (Supplementary Table 14).When comparing with Col-0 the RTW of pip1;3 and the pip2;6 RTW was higher.
In the above experiments, we studied mutants of genes with high expression levels in the roots and with suspected involvement in carbon uptake (aquaporins, sultr1;2), fixation (bca4, pepc3), and toleration of bicarbonate treatment, as well as their role in plant homeostasis (slah3).The leaf-specific aquaporin gene pip2;6 was also studied for comparison.The considerably inferior growth and photosynthetic parameters of a given mutant at 2 mM pH 5.6 NaHCO 3 compared to Col-0 could indicate a defective carbon uptake or fixation capacity, and thus, the role of the mutated gene in the process.As mentioned above, many of our mutants (sultr1;2, slah3, bca4, and pepc3) did not show a characteristic increase in ROW and RTW or photosynthetic activity as a result of bicarbonate treatment, as was expected.However, to our surprise, the pip1;3 and the pip2;6 aquaporin mutants were rather superior to Col-0 in many parameters, especially growth indicators, showing that these genes might have a negative role in the carbon uptake.Previous research e.g (Rogers et al., 1992) has shown that inorganic carbon treatment results in increased root growth and is associated with the root ultrastructure involved in carbon uptake.We validated this hypothesis and compared the growth of Col-0 and the KO mutants (Figures 7, 8; Supplementary Table 16).
Root diameter (RTD) significantly increased in Col-0 and all mutants, except for bca4, by more than 10%.The pepc3 mutant had the thinnest roots, whereas the bca4 mutant showed the smallest increase as a result of the treatment, as previously mentioned (Figure 7C).In general, the bicarbonate treatment slightly increased the epidermis width (EpWd) in Col and the mutants, with the exception of pip1;2-2, which exhibited a greater increase of approximately 46%.In the case of bca4, pepc3, and slah3 it decreased, but the difference was not statistically significant (Figure 7D).The endodermis width (EnWd) showed significant positive differences between the control and treated plants for Col-0 bca4, and slah3, with values higher than 20%.No significant differences were between the control and the treated pepc3 and aquaporin mutants (Figure 7F).The increase in the width of the cortex (CxWd) (Figure 7E) in Col-0 between the control and treatment was significant (27%).Similarly, the mutant genes pepc3 and pip1;2-2 showed a significant increase between the control and treated plants, approximately by 20%.However, no significant differences were observed among bca4, sultr1;2, pip1;2-1, and slah3 mutants (Figure 7F).The stele diameter (StDm) increased significantly in Col-0 by 14%, whereas the sultr1;2 mutant showed a substantial increase of 39%, and the pip1;2-2 mutant showed an increase of 35%.However, no significant differences were observed between the other mutants and the control and treated plants (Figure 7G).About the outer cortical cell width (OuCCWd), Col-0 showed a significant increase of 31%) between the control and the treated plants.The pepc3 and pip1;2-2 mutants significantly increased by 39% and 18%, respectively.The bca4, pip1;2-1, slah3, and sultr1.2mutants did not show significant differences (Figure 7H).Therefore, in general, many of our mutants (especially bca4, sultr1;2, slah3) showed some deviations from Col-0 in the root ultrastructure caused by bicarbonate treatment (especially in StDm, CxWd, and OuCCW), indicating that the mutated genes were missing for the uptake/fixation/ toleration of bicarbonate.

Discussion
Theoretically, root-based inorganic carbon uptake could boost photosynthetic activity by improving carbon supply.In this study, we characterized the physiology and mechanism of root-based TABLE 3 Gene candidates from the transcriptome sequencing of Arabidopsis thaliana when treated with the growth promoting system (2mM NaHCO 3 ) that might take part in root-based uptake or fixation of inorganic carbon.

Gene
Abbreviation Gene expression values and fold change (FC) of selected genes are shown in treated roots (MSRC vs. MSRT) and leaves (MSLC vs. MSLT).n.s., non-significant difference.
Gamarra Reinoso et al. 10.3389/fpls.2024.1448432Frontiers in Plant Science frontiersin.orginorganic carbon uptake by Arabidopsis thaliana (Col-0).First, a hydroponics-based inorganic carbon treatment was optimized, which served as the basis for further experiments.Inorganic carbon uptake and fixation from solutions has been proved in beans and sunflowers only at low pH (mostly in carbon dioxide form) (Shimono et al., 2019).Arabidopsis is capable of rootbased inorganic carbon uptake at a low (5.9) (Peŕez-Martıń et al., 2021) and higher pH (pH 7, 8.3) (Dabrowska-Bronk et al., 2016;Peŕez-Martıń et al., 2021).However, a higher pH might decreases the availability of some macronutrients such as phosphorus, calcium, and magnesium, and considerably hampers the uptake of micronutrients such as zinc, copper, manganese, and especially iron (Fimbres-Acedo et al., 2023).In our study, the growth promotion was achieved by supplying 2 mM NaHCO 3 at a low pH (5.6).Nutrient deficiency could explain the lack of growth promotion in the presence of MS medium and 2 mM NaHCO 3 at high pH (7.2); however, we do not exclude the possibility that high-pH growth promotion could work with different media.Another fact also argues that bicarbonate uptake could be more effective at low pH, plant roots tend to take up anions with proton cotransporters and so there is a proton requirement, though part of bicarbonate could turn into CO 2 at low pH which decreases the efficiency of the mechanism.Our candidates for bicarbonate uptake are SULTR and NRT pumps which are also known to be protoncoupled symporters (Sun and Zheng, 2015;Wang et al., 2021); and literature evidence also shows higher sulfate ion uptake with SULTRs at lower pH (Cao et al., 2013;Wang et al., 2021).
the cortex, and the role of root hairs in bicarbonate uptake (Rogers et al., Dabrowska-Bronk et al., 2016;Wang et al., 2020).Our results confirmed previous reports that the root cortex widens in response to elevated carbon supply, however, we did not observe a significant change in root hairiness.However, there was an increase in the number of lateral root branches, which could reflect a higher demand for plant growth and mineral uptake.
The beneficial bicarbonate treatment used in this study (2 mM NaHCO 3 at pH 5.6) not only improved the growth parameters but also the photosynthetic ones.In Col-0, P n was higher as a result of the carbon treatment, similar to that reported by Wang et al. (2020); however, in our study, the transpiration rate and stomatal conductance increased but decreased in Wang et al. (2020).Notably, Wang et al. (2020) used 720 ppm CO 2 different from ours since we used a soluble form of inorganic carbon added to the roots.Li et al. (2023) treated Coix lacryma-jobi with different concentrations of NaHCO 3 .Their 2 mM NaHCO 3 treatment increased the plant biomass, P n , gs, and E, similarly to our results.
Although we successfully established a hydroponics-based growth-promoting inorganic carbon treatment, it was essential to conclusively prove that the plants took up and the supplied carbon and, if possible, the ratio and thus the importance of this carbon.Carbon uptake and fixation have been studied using isotope labeling and fluxomic measurements.A higher level of 13 C aspartic acid in the xylem sap than in the phloem can be a result of rootbased PEPC-based fixation.PEPC is known to be less discriminative against the 13 C isotope compared to the green-tissue specific RuBisCo (Tcherkez et al., 2011).Miszalski et al. (2023) reported lower discrimination for 13 C in roots compared to the lamina or leaf bundle of Plantago lanceolata L. In our experiments, the concentration and ratio of 13 C-aspartic acid increased slightly when the roots were provided with 13 C. Fixation of bicarbonate by PEPC in roots has also been demonstrated by others by detecting increased oxaloacetate-derived carbonic acid (malate and citrate) and amino acid (aspartate, asparagine, lysine, methionine, threonine, and isoleucine) in the xylem sap of bicarbonate-treated roots (Bialczyk andLechowski, 1992, 1995;Yang et al., 1994;Alhendawi et al., 1997;Wanek and Popp, 2000;Nomura et al., 2006).Hibberd and Quick (2002) found evidence of C 4 -like assimilation of root-derived carbon in the stems of tobacco plants.However, the slight increase in the labeled product in our study and the decrease in PEPC gene expression (see above) led us to conclude that PEPC-based fixation in the roots could not have been the main fixation pathway of hydroponically added carbon in A. thaliana, and it would have placed huge energy demand on the roots.
On other hand, the concentration and ratio of 13 C-labeled sucrose greatly increased in the phloem sap as a result of the treatment in our work, so the transport of inorganic carbon to the shoots and Calvin-cycle-based fixation might have been of pivotal importance.Higher stomatal conductance of the treated plants could have also influenced the results as it leads to higher diffusion of 13 CO 2 , less discrimination against 13 C, and its fixation from the air through stomata (S ́liwa et al., 2019;Miszalski et al., 2023).Re-fixation of carbon originally fixed in the roots into organic molecules cannot be excluded either.Therefore, it is difficult to determine the contribution of this root-supplied inorganic carbon to the production of photoassimilates.The percentage of 13 C-labeled sucrose grew by 12-fold and reached approximately 2% of all sucrose; however, the contribution of hydroponically added 12 C could have been higher because 13 C CO 2 is known to diffuse slower and is discriminated against RuBisCo and other Calvin-cycle enzymes (Tcherkez et al., 2011).
Transcriptome sequencing was used to identify different pathways that show upregulation and downregulation during bicarbonate treatment.The expression of many more genes changed significantly in the roots (1472 DEGs) than in the leaves (341 DEGs), this may indicate the prevalence of root-based factors of inorganic carbon uptake when a source of carbon is added directly to the roots.Notably, some pathways were upregulated, as observed in earlier studies that utilized high-concentration inhibitory bicarbonate treatments.For instance, iron uptake mechanisms were activated; in fact, the ferric reduction oxidase 5 gene (AT5G23990) showed the highest fold change (FC) of all DEGs (more than 3000-fold) and the GO term "iron ion transmembrane transporter activity" and "iron ion binding" were upregulated in the roots, but "iron ion homeostasis" decreased in the leaves.High-concentration bicarbonate treatment is known to cause iron deficiency (Lucena, 2000) and upregulate iron-uptake mechanisms (Msilini et al., 2009;Peŕez-Martıń et al., 2021).Our plants might have struggled to take up iron, but with this 2 mM NaHCO 3 treatment, they were successful; they were healthy and larger than the controls.Other than FRO5, another bicarbonateprotective gene, SLAH3 (AT5G24030), was also induced in our plants.Slow-type anion channels are efflux-type nitrate and chloride channels, carbon sensors in stomata, and are involved in their closure as a result of high inorganic carbon levels (Li et al., 2022).However, SLAH3 is also highly expressed in the root xylem pericycle and pumps nitrate into the xylem, which has been shown to be a protective mechanism against high bicarbonate-induced toxicity, given the signaling role of nitrate (Duan et al., 2018).The SLAH3 mutants in our bicarbonate treatments showed somewhat reduced growth and different root structures compared with the wild type.
The low concentration of bicarbonate applied in this study and other studies may have caused beneficial and stimulating stress.In our work, carbon fixation was proven by fluxomics, and literature also suggests increased N acquisition like in tomato (Navarro et al., 2000), poplar and elder (Wanek and Popp, 2000) therefore, we believe this eustress is only partly responsible for the observed beneficial growth-promoting effect of bicarbonate treatment.The upregulation of the GO term "response to nitrogen compound" was observed in both the roots and leaves of the plant.This upregulation occurs when plants experience high levels of CO 2 , leading to an imbalance between carbon and nitrogen.Plant productivity relies on its ability to take up nitrogen and coordinate the stimulation of carbohydrate production and growth by enhancing photosynthetic processes (Li et al., 2017).The upregulation of nitrate uptake and reduction is influenced by glucose and sucrose, which increases the expression of genes involved in nitrate transport (NRT) and nitrate/ nitrite reduction (NIA/NIR) families (Thompson et al., 2017).
Since the photosynthesis, N-assimilation, and S-assimilation pathways are co-regulated, changes in one of them can influence the others (Lee and Kang, 2005).At higher CO 2 levels, upregulation of sulfate uptake and reduction was expected as happened in our study where the GO term "sulfur compound metabolic process" was downregulated but the KEGG term "cysteine and methionine metabolism" was upregulated.This predicted upregulation is a result of increased protein synthesis due to higher CO 2 levels, which leads to a greater demand for sulfur-containing amino acids such as cysteine and methionine.Cysteine is not only an important building block for proteins but can also be used to produce methionine.Methionine, similar to cysteine, is important for protein synthesis and playing a role in initiating mRNA translation.Additionally, it serves as a precursor for S-adenosyl methionine (SAM), which regulates crucial cellular processes such as cell wall biosynthesis and cell division (Hacham et al., 2013).
GO terms related to homeostasis in response to high levels of inorganic carbon in the roots and leaves suggest that the addition of bicarbonate to the root environment could have affected internal pH and nutrient availability, potentially causing osmotic, ionic, and oxidative stress (Li and Yang, 2023).To sustain plant growth, ion and pH homeostasis must be maintained.The upregulation of GO terms, such as the 'response to oxidative stress', 'cation and ion homeostasis', 'cellular homeostasis', and 'chemical homeostasis' in the roots, reflects the presence of these stresses.Cellular pH is tightly regulated to ensure optimal conditions for biochemical reactions, and plant cells possess numerous transport systems comprising cation/H + and anion/H + exchangers that maintain pH homeostasis in different subcellular compartments as well as transmembrane H + gradients as a proton motive energy source (Rottet et al., 2021).
The GO term 'trehalose biosynthetic process' was upregulated both in the roots and leaves.The trehalose-6-phosphate (T6P) signaling pathway regulates the relationship between sources and sinks under different environmental conditions.This regulation leads to an increase in soluble saccharide content, resulting in a significant increase in aboveground weight, plant height, and stem diameter of the plant (Li et al., 2024).T6P has emerged as an important signaling metabolite that regulates carbon assimilation and sugar status in plants, plays an essential role in plant development, and is an important structural component of plant cells (Ponnu et al., 2011).Changes in photosynthetic capacity in plants with altered T6P levels are thought to be influenced by T6P's impact on leaf development, particularly on cell division and cell wall biosynthesis.In addition to its role in regulating carbohydrate metabolism, T6P is also important for normal plant (Ponnu et al., 2011).T6P is also believed to play a role in upregulating multiple macronutrient pathways when carbon levels are high, as demonstrated by Wingler et al. (2012).
Interestingly, the expression of high-level b-carbonic anhydrases (bCA1 and bCA4) did not increase as a result of bicarbonate treatment in the transcriptomic analysis.The level of bCA4 significantly decreased in the roots, and the expression was higher in the leaves.This contrasts with the results of Dabrowska-Bronk et al. (2016), who found increased expression of many bCA genes, especially bCA1 and bCA4, following bicarbonate treatment, although the experimental conditions were different from our research (hydroponically grown plants at pH 5.6 vs. irrigating potted plants with a pH 7 solution).However, the bca4 mutants clearly showed growth differences from the wild type in our bicarbonate treatments; therefore, we agree with the authors that "CAs allow to use of soil HCO 3 − ions as a substrate for photosynthesis in foliar tissues" (Dabrowska-Bronk et al., 2016).
From the various mechanisms, that we hypothesized would participate in root-based carbon uptake, many aquaporin channel genes were activated as a result of the treatment.Some aquaporin genes, such as the well-characterized carbon channel AtPIP1;2 (Heckwolf et al., 2011) excelled with very high expression levels but were only slightly induced further because of the treatment (FC = 1.33).On the other hand, many nitrate, phosphate, and sulfate anion channels, both high-and middle-level channels (e.g., SULTR1;1, SLAH3, NRT2.1, and PHO1) were significantly induced (FC of 2-11).It has been hypothesized that some anion transporters may also have some affinity for bicarbonate [reviewed by Poschenrieder et al., 2018].Therefore, increased expression of these anion transporters might take up bicarbonate ions, or they might be required for the uptake of other nutrients to sustain increased growth.As was already mentioned, we hypothesize that SULTRs and NRTs take up bicarbonate ions into the root, while PHO1 as an efflux-type phosphate channel pumps, phosphate into the xylem (Stefanovic et al., 2011) and it might also pump bicarbonate.
Among anion channels, sulfate transporters have been implicated in having some affinity for bicarbonate (Frachisse et al., 1999;Poschenrieder et al., 2018;Wang et al., 2021).Arabidopsis has two high-affinity sulfate transporters (SULTR1;1 and SULTR1;2) that represent the sulfate uptake activity at the root surface (Maruyama-Nakashita et al., 2003).SULTR1;2 is predominantly expressed in the root cortex, involved in importing sulfate from the environment into the root, and is not a very specific transporter; it transports selenate and chromate as well (Shibagaki et al., 2002;Maruyama-Nakashita et al., 2003;Xu et al., 2021).In our experiments, the high expression level of the Proposed root-based inorganic carbon uptake and fixation mechanisms of Arabidopsis thaliana.Bicarbonate (HCO 3 -) ion is taken up by anionproton symporters, a minority of bicarbonate is fixed in the roots by PEPC (phosphoenulpyruvate carboxylase), but the majority is transported up to the shoots by the xylem and fixed there in the Calvin-cycle after conversion to CO 2 by carbonic anhydrases (CAs).Part of bicarbonate may be converted to CO 2 by CAs and vented out of the root cells through aquaporin channels.However CO 2 may leak into the xylem at root tips and branching sites where the Casparian strip has not formed.Note the pH of the different spaces (root/leaf symplast: alkaline, apoplast: acidic, xylem sap: acidic) which also affects carbon fixation and transport, inorganic carbon can be found predominantly as CO 2 at lower and bicarbonate at higher pH but the conversion is slow without CAs.OAA stands for oxaloacetate, TCA, tricarboxylic acid cycle; E, epidermis; C, cortex; En, endodermis and S, stele.
sulfate transporter SULTR1;2 was slightly further induced (FC = 1.44), the lower level SULTR1;1 was greatly induced (FC = 3.57) as a result of bicarbonate treatment.Under normal conditions with adequate sulfate in the soil, SULTR1;2 is the primary transporter for sulfate uptake.However, under stressful conditions, SULTR1;1 is the primary sulfate uptake transporter (Yoshimoto et al., 2002).
To determine whether aquaporin channels or anion pumps are more likely to be responsible for inorganic carbon uptake by roots, and the role of bCA and PEPC, gene-knockout mutants were tested for growth promotion under bicarbonate treatment.The sultr1;2 mutant showed the greatest deviation (inferior growth promotion and photosynthetic activity) compared to the treated wild-type, which underlines the importance of this gene as a possible bicarbonate uptake pump.Treated roots exhibited various ultrastructural changes.While the growth and photosynthetic activity of many mutants were inferior to those of Col-0, the aquaporin mutant pip1;3 showed superior performance; therefore, we believe that this gene could play a negative role in root-based carbon uptake.Based on our findings (transcriptome and mutants), we hypothesize that despite the low pH (carbon-dioxide-rich environment), Arabidopsis roots still preferentially take up bicarbonate ions with anion pumps, with possible candidates being SULTR1;1, SULTR1;2 and NRT2;1, while PHO1 might transport the bicarbonate into the xylem.The bicarbonate ion predominantly stays in this form in the root symplast, as the pH is usually slightly alkaline in the cytosol of plant cells, ranging from 7.2 to 7.5 (Felle, 2001).Carbon dioxide-permeable aquaporins may have a negative effect by venting the carbon dioxide out of the plant cells, produced by carbonic anhydrases.We found that bCA4 expression decreased in the roots, possibly preventing the loss of CO 2 .The growth and photosynthetic activity of the pepc3 mutants could not be boosted as much by our bicarbonate treatments as in the case of the Col-0.As it was mentioned, the production of labelled aspartic acid only slightly increased in the xylem sap of 13 Ctreated plants.Large-scale production of malate, another possible derivative of PEPC-based oxaloacetate, is yet more unlikely in the roots due to the requirement of reductive power (one NADPH corresponds to the energy of three ATPs), although the production of small amounts of malate and other organic acids has been detected in carbon-treated roots (Cramer et al., 1993;Alhendawi et al., 1997).However, PEPC-based carbon fixation can play a role in N acquisition; some articles refer to better ammonia assimilation in bicarbonate treated-roots (Marıń-Peña et al., 2024), and others to nitrate, (Viktor and Cramer, 2005) in the roots, or both nitrate and ammonia (Andrews et al., 2019).Apart from the proposed bicarbonate uptake mechanism through proton-anion symporters, it can also be anticipated that CO 2 (in this form, as the apoplastic pH is acidic in general: Felle, 2001) leaks into the xylem at root tips and branching sites where the Casparian strip has not formed.The xylem sap pH is generally acidic in Arabidopsis (Grunwald et al., 2021), therefore presumably mostly CO 2 is transported up to the shoots in the xylem sap, enters the mesophyll cells, where it might be converted to bicarbonate at the alkaline cytoplastic pH but converted back to CO 2 by CAs for fixation for RuBisCo.
Our study found growth-promoting effects of low-concentration bicarbonate treatment in Arabidopsis with root-based uptake.Higher growth was possibly regulated by trehalose-6-phosphate signaling, and various mechanisms have been proposed to support the increased growth, including anatomical changes (wider root cortex) and higher nitrogen and sulfur assimilation.Transcriptomic data revealed the involvement of various transmembrane transporter activities, cation and anion transport, and maintenance of homeostatic balance.Increased iron uptake and higher SLAH3 expression alleviated the toxic effects of bicarbonate.In addition, we identified potential inorganic carbon transporters in the roots that can take up either HCO 3 -.Further studies are necessary to confirm these findings.Our investigation also revealed that unlike in leaves, aquaporins, especially PIP1;3 may have a negative impact by venting carbon dioxide out of the plant in the roots.These results suggest that a minority of bicarbonate was fixed in the roots by PEPC, but the majority was transported up to the shoots by the xylem and fixed there by RuBisCo after conversion to CO 2 by carbonic anhydrases.The proposed root-based inorganic carbon uptake and fixation mechanisms of Arabidopsis thaliana are presented on Figure 9.The study shed new light on the mechanisms of root-based inorganic carbon uptake, which could be relevant for boosting photosynthesis with carbon supply.0205K0036P]'.LG, JJ and EJ are grateful for Stipendium Hungaricum scholarship from Hungarian government.

FIGURE 1
FIGURE 1The concentration of aspartic acid and sucrose (nM) in the xylem and phloem sap of plants treated with 0mM, 12 C 2 mM NaHCO 3 and 13 C 2 mM NaHCO 3 as sources of inorganic carbon applied to the roots.(A) Level (nM) of 12 C, 13 C aspartic acid in the xylem and phloem sap, (B) Level (nM) of 12 C, 13 C sucrose in the xylem and the phloem sap.Each value represents the mean ± SD of n=6.Asterisks indicate statistically significant differences (p-value ≤0.05*, ≤0.01 **, ≤0.001 ***, Paired Student t-test).

FIGURE 2
FIGURE 2 Volcano plots and illustration of differentially expressed genes in MSLC vs MSLT and MSRC vs MSRT.(A) Representation of the volcano plot for the MSLC vs MSLT.(B) Representation of the volcano plot for the MSRC vs MSRT.(C) The total, the downregulated, and the upregulated DEGs are presented in the MSLC vs MSRT and MSRC vs MSRT.The volcano plots and the bar representation were graphic considering a based Log 2 Fold change ≥ 1/-1 and P-adj ≤ 0.05.

FIGURE 3
FIGURE 3 Venn diagrams and bar plot comparing MSLC vs. MSLT and MSRC vs. MSRT.A Venn diagram of all the DEGs present in MSLC vs MSLT and MSRC vs MSRT.(B) Venn diagram of upregulated DEGs present in the MSLC vs MSLT and MSRC vs MSRT.(C) Venn diagram of downregulated DEGs present in the MSLC vs MSLT and MSRC vs MSRT.(D) The number of DEGs specifically regulated by MSLC vs MSLT and MSRC vs MSRT.

FIGURE 4 GO
FIGURE 4 GO analysis for the leaves and the roots.(A) Upregulated GO terms for MSLC vs MSLT.(B) Downregulated GO terms for the MSLC vs MSLT.(C) Upregulated GO terms for the MSRC vs MSRT.(D) Downregulated terms for the MSRC vs MSRT.

FIGURE 7
FIGURE 7 Growth parameters and root microscopy analysis of Col-0 and the mutants in the roots treated with 0mM NaHCO 3 and 2mM NaHCO 3 .Root length (A), root weight (B), root diameter (C), epidermis width (D), Cortex width (E).(F) Endodermis width (F).Stele diameter (G).Outer cortical cell width (H).Each value represents the mean ± SD of n=12 for growth parameters and the root analysis.Asterisks indicate statistically significant differences between the control and the treated plants (p-value ≤0.05*, ≤0.01 **, ≤0.001 ***, Student t-test).

TABLE 1
Transporters and channels upregulated in the MSRC vs. MSRT.

TABLE 2
Transporters and channels downregulated in the MSRC vs. MSRT.comparison of the mutants and the Col-0 in the P n values and transpiration rates of the pip2;6 had the lowest values within all the mutants., gs.Representative images of the controls (Figures